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ABSTRACT 


Displacement equations are developed for a cantilever tubular wing spar under bending, torsion, and 
combined bending and torsion loading. The displacement equations are expressed in terms of strains 
measured at multiple sensing stations equally spaced on the surface of the wing spar. The bending and 
distortion strain data can then be input to the displacement equations to calculate slopes, deflections, and 
cross-sectional twist angles of the wing spar at the strain-sensing stations for generating the deformed 
shapes of flexible aircraft wing spars. The displacement equations have been successfully validated for 
accuracy by finite-element analysis. The strain-displacement theory that has been developed could also be 
applied to calculate the deformed shape of simple and tapered beams, plates, and tapered cantilever wing 
boxes. The displacement equations and associated strain-sensing system (such as fiber optic sensors) form 
a powerful tool for in-flight deformation monitoring of flexible wings and tails, such as those often 
employed on unmanned aerial vehicles. Ultimately, the calculated displacement data can be visually 
displayed in real time to the ground-based pilot for monitoring the deformed shape of unmanned aerial 
vehicles during flight. 


NOMENCLATURE 

A cross-sectional area of uniform cantilever beam, in 2 

a mean radius of four-ply composite tube wall and aluminum tube wall, or mean radius of Helios 
31 -ply reinforced spar cap wall, in. 

c depth factor of beam (distance between neutral axis to bottom surface of beam), in. 
deg degree 

E Young’s modulus, lb/in 2 

Ei Young’s modulus of lamina in fiber direction, lb/in 2 

Ej Young’s modulus of lamina in direction transverse to fiber direction, lb/in 2 

e finite-element span-wise length, in. 

G shear modulus, lb/in 2 

Ci m shear modulus of lamina, lb/in 2 

I moment of inertia of tube cross section, in 4 
i identification number for strain-sensing station (=0, 1, 2, 3,. . ., n ) 

J polar moment of inertia, in 4 

j index 

/ length of cantilever tube, in. 

M moment, in- lb 

Mj moment at strain-sensing station, X ( , in- lb 
n identification number for the last strain-sensing station, x n (= /) 

P applied load, lb 

SPAR Structural Performance And Resizing 

Tj twisting moment at distortion-sensing station, Xj , in-lb 

t thickness of four-ply composite (or aluminum) wall, in. 

t thickness of 3 1 -ply composite reinforced wall region, in. 



UAV unmanned aerial vehicle 
x, y Cartesian coordinates, in. 

X[ i -th strain-sensing station, or its axial coordinate, x = X [ , in. 
y deflection in y-direction at any point, x, between two adjacent strain-sensing 
stations, Xj_ \ < x < Xj , in. 

yi deflection in y-direction at strain-sensing station, Xj , in. 

Yi surface distortion angle (shear strain) at distortion-sensing station, Xj , rad 
A l = l/n , distance between two adjacent strain-sensing stations (equally spaced), in. 

(5] , 82 displacement at two ends of finite element, in. 

£ bending strain at any point, in/in. 

£[ bending strain at strain- sensing station, Xj , due to bending only, in/in. 

£[ true bending strain at strain-sensing station, jc; , due to bending and torsion, in/in. 

£■ principal tensile strain at distortion-sensing station, x/ , in the 45-deg helical direction, in/in. 

9 beam slope at any point, x, between two adjacent strain-sensing stations, x/q < x < x,- , rad or deg 

Oj beam slope at strain-sensing station, x, , rad or deg 

V Poisson’s ratio 

V LT Poisson’s ratio of lamina 

a bending stress of outermost fiber, lb/in 2 

<T- principal tensile stress at twist-sensing station, x/ , in the 45-deg helical direction, in/in. 

(Tmax )/ maximum shear stress at distortion-sensing station, i, lb/in 2 
(f>l cross-sectional twist angle at strain-sensing station, X [ , rad or deg 
( )i+ outboard side of strain-sensing station, x,- 
( )/_ inboard side of strain-sensing station, Xj 


INTRODUCTION 


Many unmanned aerial vehicles (UAVs) currently in service today must be ultralightweight and 
designed with highly flexible structural control surfaces. One such vehicle, the Helios flying wing, has a 
wingspan of 247 ft. The wingspan of the Helios is longer than the wingspans of the U.S. Air Force C-5 
military transport (Lockheed Martin Aeronautics Company, Marietta, Georgia), which is 222 ft, and the 
Boeing 747 commercial jetliner (Boeing Company, Seattle, Washington), which ranges from 195 to 215 
ft, depending on the model. Because of an extremely long wingspan, the wings of vehicles such as the 
Helios will undergo extra-large deformations into high dihedrals during flight, with wingtip deflections 
reaching as high as 40 ft. This type of deformed shape is similar to a 72-deg circular arc with a 197-ft 
radius. 

Recently the Helios wing broke up in midair while flying at an altitude of approximately 3000 ft. The 
primary cause of the mishap is believed to have resulted from undamped pitch oscillations of the highly 
deformed wing, resulting in the breakup of the wing. To avoid future mishaps, a method must be 
developed in which the wing deformation can be visually monitored during flight. This way, if the 


2 



wingtip deflections approach the design limit, the ground-based pilot could execute timely emergency 
maneuvers to avoid further deformations. 

The major load-carrying main spar of the Helios wing is a carbon fiber composite tube (4-in. radius), 
which is thicker on the top and bottom (spar caps) to improve bending stiffness. The spar also is wrapped 
with Nomex" and Kevlar® (both registered trademarks of E.I. duPont de Nemours and Company, 
Wilmington, Delaware) for additional strength. If multiple strain sensors, such as lightweight fiber optic 
strain-sensing systems, are installed at discrete sensing stations on a beam-like structure (such as the long- 
span highly flexible Helios flying wing), the strain sensor data can be used to calculate the deflections and 
cross-sectional twist angles of the beam-like structure and thereby predict its deformed shape during 
flight. 

This report discusses the development of theoretical displacement equations for uniform cantilever 
beams subjected to bending, torsion, and combined bending and torsion. The displacement equations are 
expressed in terms of strains, the output of multiple strain sensors embedded on the wing spar surface, for 
predicting wing deflections and cross-sectional twist angles during flight. The displacement equations and 
onboard strain-sensing system could form a powerful tool for in-flight wing deformations monitoring by 
the ground-based pilot for maintaining safe flights of highly flexible flying wings. 

HELIOS UNMANNED AERIAL VEHICLE (UAV) 

FLEXIBLE STRUCTURE 

The Helios prototype (fig. 1) is an ultralightweight, unmanned, solar-powered flying wing aircraft 
(1,600-lb weight) designed to fly at altitudes up to 100,000 ft. Fourteen electric motors (2 hp each) are 
required to drive 79-in. diameter composite propellers and other systems. The power required for these 
motors is generated by high-efficiency solar cells spread over the upper surface of the wing. The cruising 
speed of the Helios is in the range of 19 to 27 mi/h. 

The Helios wing has an 8-ft chord length uniform over the 247-ft wingspan. As discussed previously, 
the wingspan of the Helios is longer than the wingspans of the U.S. Air Force C-5 military transport (222 
ft) and the Boeing 747 commercial jetliner (195-215 ft, depending on the model). The maximum wing 
loading on the Helios 1,976-ft 2 wing is 0.81 lb/ft 2 . 

IN-FLIGHT WING SHAPE MONITORING 


In-flight wing deformation monitoring of ultralightweight flying wing aircraft such as the Helios is 
particularly difficult because of weight restrictions and highly flexible nature of the extremely long-span 
structure. Several methods for in-flight monitoring of wing deflections exist for use on conventional 
aircraft. 

One method is the electro-optical flight deflection measurement system, which is composed of 
onboard optical receivers and several wing-mounted light- emitting targets (ref. 1). This system provides 
wing displacement information during flight but is too heavy for lightweight flying wing applications. 

Another in-flight deflection measurement method is the use of conventional strain gages to collect 
local strains for wing deflection information. Numerous strain-sensing stations are required to capture the 
higher displacement modes of the flexible wings. When multiple strain-sensing stations are used, the 
weight of the strain gage lead wires alone is too heavy and impractical for most weight-conscious 
lightweight flying wing aircraft. 
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Fiber optic sensors offer the most attractive alternative to conventional strain gages. They are 
lightweight, fine, and flexible filaments (approximately the size of human hairs), and they can be highly 
multiplexed at desired sensing intervals. 

THEORETICAL APPROACH 


If multiple strain sensors are installed at discrete sensing stations on a cantilever wing spar, sensor 
data can be used to calculate the deflections and twists of the wing spar during flight. Classical beam 
theory can be used to develop theoretical slope, deflection, and cross-sectional twist angle equations for 
the uniform cantilever beam. These displacement equations are written in terms of the strains for 
calculating the deformed shape of the beam. In the present report, the input strains are calculated from the 
displacement output of finite-element stress analysis. The accuracy of the newly developed displacement 
equations (for bending and torsion) for the uniform cantilever beam is then verified by the finite-element 
displacement output. With slight modifications, the displacement equations developed for the uniform 
cantilever beam could be applied to predict the deformed shape of simple beams (two-point supported), 
tapered beams, plates, and tapered cantilever wing boxes. 

UNIFORM CANTILEVER BEAMS 


In this section, basic displacement equations are developed for a uniform cantilever beam (for 
example, the tubular spar of the Helios wing) subjected to bending, torsion, and combined bending and 
torsion. Although the highly flexible wing of the Helios is capable of undergoing extraordinarily large 
deformations (wingtip deflections up to 40 ft), the local strains on the wing spar surface remain small. 
Therefore, the classical small strain theory can be used to derive the slope and displacement equations for 
the uniform cantilever beam. 


Bending 

In this section, displacement equations for bending only are developed. These equations are used to 
describe the local slopes and deflections at points along the strain-sensing line on the bottom surface of a 
uniform cantilever beam of arbitrary cross-sectional shape (for example, circular tubes of the Helios spar). 

Moment-Strain Relationship 

The classical bending equation for the uniform beam is given by (refs. 2, 3) 

d 2 y _ M(x) m 

dx 2 El 

in which y is the vertical displacement, x is the span-wise coordinate, M(x) is the bending moment, E is 
the Young’s modulus, and / is the moment of inertia. 
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The bending moment, M(x ) , in equation (1) may be related to the associated bending strain, e(x ) , at 
the bottom or top fiber (for example, generatrix of the Helios spar tube) through the outermost fiber 
bending stress equation, o(x) = [M(x)c\/l (refs. 2, 3): 


£(x) = 


<r(x) 

E 


M(x)c 

El 


a(x) = 


M(x)c 

I 


( 2 ) 


In equation (2), c is the depth factor of the beam (distance between the neutral axis and the bottom 
surface of the beam). In light of equation (2), the beam bending equation (1) can be written in term of 
£(x) as 


d 2 y _ £(x) 

j 2 l -3 ' 

ax c 

Note that under the strain formulation, beam bending equation (3) contains only the beam geometrical 
parameter, depth factor, c, and that the flexural rigidity term, El, is eliminated. 

Consider the case of a highly flexible, extremely long tubular spar of the Helios wing that carries 
concentrated weight (solar-powered motors, pods, hydrogen storage tanks, and so forth) at different wing 
stations. Because the wing is subjected to aerodynamic lift forces during flight, the wing spar tube 
bending moment, M (x ) , could be a complicated nonlinear function of the span-wise coordinate, x. 

Strain-Sensing Stations 

Figure 2 shows the uniform cantilever tubular spar with a length of l and an outer radius of c. The 
beam is installed with n + 1 equally spaced bending strain sensors at the sensing station, Xj (i = 0, 1, 2, 3, 
..., n), along the bottom or top generatrix, which is a straight line on the tubular surface parallel to the 
axis, x, of the tubular spar. The first and final strain-sensing stations { xq (x=0) , x n (x-l ) } are located at 
the built-in end and the beam tip, respectively. The beam is now discretized into n sections with equal 
length, Al = Un . Although the bending strain is zero at the spar tip, the «th bending strain sensor is 
installed there because of mathematical convenience in the derivation of the displacement equations. 

The bending moment, M(x), can be assumed to be a piecewise-linear function along the axial 
coordinate, x. Namely, in the region, x/_] < x < x/ , between any two adjacent strain-sensing stations, 
{ X;_i , Xj }, M(x) is considered as a linear function of (x- x/_] ) as 

M(x) = M i _i-(M i _i~M i ) X X '~ l ; x M <x<x; (4) 

Al 

in which are the bending moments at the two adjacent strain-sensing stations, (x/_| ,x ; - } , 

respectively, and M ( = x ; - — x;_i = l/n ) is the axial distance between the two adjacent strain-sensing 
stations, (x/_| ,x; } . 
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In light of equation (2), local moment equation (4) can be written in terms of the local bending strain, 
e(x) , for the region, x- t _\ <x< Xj , as 


X ] 1 ; Xi_i<x<Xj ; Al = - ( 5 ) 

Al n 

in which { £/_] , £,• } are the bending strains measured at the two adjacent sensing stations, { x,-_ | , x/ } , 
respectively. 

Slope Equations 

The slope, tan 0(x) , of the uniform beam in the region, x/_] < x < x/ , between the two adjacent 
strain-sensing stations, { x/_ | , x/ ) , can be obtained by integrating equation (3), with the constant of 
integration determined by enforcing the continuity of the slope at the adjacent inboard strain-sensing 
station, x/_ ] , as 


,2 

rx d y rx 

tan0(x) = — -r-dx + tan0,_i = 

J X; i j2 v 1 X; 


£(x) 


dx + tan 0, 


Jx i - 1 dx 
Slope increment 


Slope at x ( -_ | 


h'-l 


'i-1 


X;. 


V_i < x < x t 


( 6 ) 


in which tan0 ; _j (constant of integration) is the slope at the adjacent inboard strain-sensing station, 
X;_ i . Substituting strain equation (5) into slope equation (6), and carrying out the integration, one obtains 
the slope, tan 0/ = [tan Q(x - t )] , at the strain-sensing station, x; , as (ref. 4) 


tan Oj = 


A l_ 
2c 


(£/_l + £j ) + tan 6[_ i 


(/' = 1, 2, 3, ..., n) 


(7) 


Equation (7) is a recursion formula showing that the slope, tan 0 ; - , at the strain-sensing station, x ; - , is 
related to the slope, tan 0 ; _ ] , at the adjacent inboard strain-sensing station, x ( -_] . Applying the 
descending indices relationship causes slope equation (7) to become 


* Al 
tan 0; = — 

2c 


i — 1 

£ 0 + 2 X £ j + £ i 
7=1 


(/ = 1, 2, 3, ..., n ) 


( 8 ) 


Note that the slope, tan 0/ , is determined from the geometrical and strain data evaluated at all the 
inboard and current strain-sensing stations, (xq , X] , xj , •••, x/ ). The first objective is to calculate tan Oj (i 
= 1, 2, 3, ..., n ) for all strain-sensing stations. If the last strain-sensing station, x n , is located at the 
wingtip, then the slope, tan 0 n , becomes the wingtip slope. 

Deflection Equations 

The deflection, y(x), of the uniform beam in the region, X;_j < x < x/ , between the two adjacent 
strain-sensing stations, { x,_j,x ; - }, can be obtained by integrating slope equation (6) with the constant of 
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integration determined by enforcing the continuity of deflection at the adjacent inboard strain-sensing 
station, x,_j , as 


rx 

v(x) = tan 6(x)dx + 

v 

Integration of slope 


-VM 

Deflection 
at Xj_i 


rx rx £(x) , , rx 

! dxdx+ tant i;_\dx + 

J Xi-\ C J Xj - 1 

V V 

Deflection increment Deflection at x 

due to tan 0j_ \ 


y^A 

Deflection 
at Xj_i 


Xj_l<X<Xi 


(9) 


in which y ; -_j (integration constant) is the deflection at the adjacent inboard strain-sensing station, x/_] . 
Carrying out the integration of equation (9) in light of strain equation (5), one obtains the deflection, 
y'i = [y(Xj )] , at the strain-sensing station, Xj , as (ref. 4) 

2 

Vi = ( ^—( 2£ i-l +£ i) + yi-l +Altane i-l > 0'= h 2, 3,..., n) (10) 

6 c 

Deflection equation (10) is also a recursion formula showing that the deflection, y \ , at the strain- 
sensing station, X [ , is expressed in terms of the deflection, y;_] , and the slope, tan , at the adjacent 
inboard strain-sensing station, x/_ | . Applying the descending recursion relationship for v/_ | in light of 
equation (10), and applying expression (8) for tan0 ; -_i , one could express deflection equation (10) for >7 
as a summation of all the inboard strains including the strain at the current sensing station, X [ : 


>V = 


(AI ) 2 
6 c 


i — 1 

(3/ - 1)£q +6^0'- j)£ j + £/ 
7=1 


(/ = 1, 2, 3,..., «) 


( 11 ) 


By plotting the deflections, y[ (i = 1, 2, 3,..., n), at the strain-sensing stations, x,-, one can construct the 
deformed shape of the tubular wing spar, which could be displayed on the screen so that the ground-based 
pilot can view the real-time wing deformation status. 

The contribution of deflection caused by shear effect usually is very small for a long cantilever beam. 
If shear effect is desired, however, an additional shear effect term must be added to equation (11) (see 
Appendix B of reference 5). 


Torsion 

In this section, displacement equations are developed for a tubular cantilever spar under torsion. 
Consider a cantilever tubular spar of length /, with n number of equally spaced distortion-sensing stations 
(fig. 2), subjected to span-wise varying twisting moments. 
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Torque-Twist Relationship 

Let Tj be the twisting moment at the distortion-sensing station, X ( , and let y ; - be the associated 
surface distortion angle (shear strain) at the distortion-sensing station, x ; - (fig. 3). The classical torque- 
distortion 

relationship can then be written as (refs.l, 2) 


7 i = 


Ti c_ 
GJ 


( 12 ) 


in which G is the shear modulus, J is the polar moment of inertia, and c is the outer radius of the wing 
spar. The surface distortion angle, y ; - (shear strain), in equation (12) is obtained from the distortion 
sensor output using 


= (W)» = °J- = l e p = 2(1 + v)ef (13) 

1 G G G 1 1 

in which (T max ) ; - is the maximum shear stress at the distortion-sensing station, i, and <7- [=(T max ); ] 

p 

and £j are the principal tensile stress and strain, respectively, in the 45-deg helical direction at the 
distortion-sensing station, i (fig. 3). 

Cross-Sectional Twists 

The total cross-sectional twist angle, (pj , at the strain-sensing station, x ; - , can be obtained by 
summing up all surface distortion angles, y j , at all the inboard strain-sensing stations, j (= 1, 2, 3,..., i - 
1), as (fig. 3) 


<Pi = — (7o + r, + 72 + 73 + + 7- 1 ) (14) 

c 

Notice that when equation (14) is written for the beam tip (/ = n), the resulting expression of the beam 
tip cross-sectional twist angle, (j ) n , does not contain the tip distortion term, y n . Therefore, installation of 
the tip distortion sensor for measuring y n is not required (fig. 2). 

Combined Bending and Torsion 

In this section, displacement equations are developed for a cantilever tubular spar under combined 
bending and torsion. The calculations of true bending strains are presented, from which modified slope 
and deflection equations are established. 

True Bending Strain 

As shown in figure 3, if the tubular spar is subjected to torsion in addition to bending, the bending 
strain-sensing station, x ; - (for measuring bending strain, £/ ), will move tangentially by a cross-sectional 
twist angle, (pj , and the bending strain-sensing axis will be tilted by a local surface distortion angle, y ; - . 
Thus, the strain output, £,■ , obtained at the sensing station, x; , no longer gives the true bending strain. To 
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obtain the true bending strain, £,• , the value of the measured strain, £/ , must be corrected by using the 
bending strain correction equation, 


cos (j)j cos Ji 


(15) 


in which the surface distortion angle, Yi ■> ar) d the cross-sectional twist angle, 0/ , are calculated from 
equations (13) and (14), respectively. 


Modified Slope and Deflection Equations 

With the effect of torsion, slope equation (8) and deflection equation (11) for the bending case must 
be modified by replacing £,• with £; : 


tan dj = 


Al 

2c 


i—1 

^0 + 2 X £j + £i 
7=1 


(16) 


>V = 


(Al) 2 

6c 


i — 1 

(3i — 1 )£q + 6^ (i-j)£j +£,• 
7=1 


(17) 


When the true bending strain data, (£q,£i, £ 2 , £3 , are substituted, the slopes and deflections of the 

tubular spar at the strain-sensing stations can be calculated from equations (16) and equations (17), 
respectively, for the combined bending and torsion loading case. 


FINITE-ELEMENT ANALYSIS 


To test the prediction accuracy of deflection equation (11), strain sensor data and measured deflection 
data are needed. Before the experimental data are available, an alternative approach is to use finite- 
element analysis to generate the desired bending strain and deflection data. The Structural Performance 
And Resizing (SPAR) finite-element computer program (ref. 6) was used for this purpose. The following 
section describes the method for generating the bending strains, £[ , from the SPAR output. 

Generation of Bending Strains 

From the SPAR nodal displacement output, the axial displacements, (5] , 82 } , at the two ends of the 
element lying in the lower outermost surface can be used to calculate the bending strain for that particular 
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element. If the strain-sensing station, X ( , is located in the middle of the element, the displacement 
differential, (52 -<5i), can be divided by the element axial length, e, to generate the desired bending 
strain, £j , at the sensing cross section, Xj . Namely, 


^2-V 

V e ); 


(18) 


If the sensing station, jc;, is located at the juncture of two adjacent elements, then the following 
average strain value is used for the strain, £j , at the sensing station, X ( : 


y$2-«ii 

+ 

r«2-«ij 


\ e ) 

i~\~ 

{ e ) 

i— _ 


(19) 


In equation (19), {()/+, ();_} denote the elements on the outboard and inboard sides, respectively, of 
the strain-sensing cross section, x ; - . 

The SPAR bending strains calculated from equations (18) and (19) can be input to deflection equation 
(11) to calculate the beam deflections. These calculated deflections can then be compared with the 
deflections obtained from the SPAR displacement output to check the prediction accuracy of deflection 
equation (11). 


Cantilever Tube Models 

Figure 4 shows a typical finite-element model for a cantilever tube with a length of l =100.5 in. The 
model consists of 100 four-node shell elements in the axial direction and 36 in the circumferential 
direction. The end disk was added so that the tip load can be applied at the disk center. The cantilever 
tube is fixed in space at the left end and is subjected to an upward vertical load of P = 100 lb at the tip. 
The following three types of tubes, fabricated with different materials, were considered: 

1 . Aluminum tube 

2. Four-ply composite tube: carbon fiber composite wall 

3. Helios composite tube: carbon fiber composite wall with additional carbon fiber composite 
reinforcements on the top and bottom of the spar caps. 

Table 1 lists the dimensions and table 2 lists the material properties of these three types of tubes. 

Table 1. Dimensions of three types of tubular spars. 


Tube type 

/, in. 

a, in. 

c, in. 

t, in. 

Aluminum 

100.5 

4.0 

4.01148 

0.02296 

Four-ply composite 

100.5 

4.0 

4.01148 

0.02296 

Helios composite 

100.5 

4.0* (4.07749)** 

4.16646** 

0.02296* (0.17794)** 


* Four-ply regions ** Reinforced regions 


10 



Table 2. Material properties of carbon fiber composite and aluminum tubes. 


Carbon fiber composite 

Aluminum 

E l = 16.3 xlO 6 lb/in 2 

E = 10.5 x 10 6 lb/in 2 

E t =0.98 xlO 6 lb /in 2 

NA 

G lt =0.71 lb/in 2 

G = 4.0 xlO 6 lb/in 2 

vjjf = 0.34 

v = 0.33 


NUMERICAL EXAMPLES - CANTILEVER TUBES 


The SPAR displacement output and strain equations (18) and (19) were used to calculate the bending 
strain, £,• , for the sensing stations in all cases. These strain values were then input to the appropriate 
deflection equations to calculate the deflections at all the strain-sensing stations of all the beam cases. 
This section presents examples of three types of cantilever tubes (aluminum, four-ply composite, and 
Helios composite). The deflections calculated from the Ko displacement theories are compared with those 
calculated from the SPAR finite-element analysis. 

Aluminum Cantilever Tube 

Figure 5 shows the bending strains, £/ , at the sensing stations on the aluminum tube, calculated from 
the SPAR nodal displacement output. The strain data from figure 5 were used to calculate the deflections, 
y'i , from deflection equation (11). Table 3 compares the deflections, >’,■ , calculated from deflection 
equation (11) with those obtained from the SPAR displacement output for the case of n = 8 (»+ 1= 9 
strain-sensing stations). Note that the difference between the deflection values calculated from the SPAR 
program and those calculated from equation (11) is minimal for the present uniform cantilever aluminum 
tube. 

Table 3. Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (1 1); aluminum cantilever tube; n = 8. 


Deflection ( y ; - ), in. 



F0 

>’l 

T2 

F3 

T4 

F5 

F6 

yi 

F8 

SPAR 

0.0000 

0.0179 

0.0636 

0.1341 

0.2253 

0.3330 

0.4530 

0.5814 

0.7138 

Df-flprllnn 

equation (11) 

0.0000 

0.0157 

0.0605 

0.1302 

0.2205 

0.3273 

0.4465 

0.5738 

0.7051 

Difference , 
percent 

0.0000 

12.2142 

4.9356 

2.9223 

2.1262 

1.7208 

1.4524 

1.3072 

1.2230 


The data in table 3 are plotted in figure 6. The two sets of deflection curves are quite close, indicating 
the accuracy of deflection equation (11). The deflection curves calculated from equation (11) lie slightly 
below those obtained from the SPAR output. 
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Four-Ply Composite Cantilever Tube 


Figure 7 shows the bending strain curve, £, , calculated from the SPAR output, for the strain-sensing 
stations on the four-ply composite tube. In this case the SPAR strain curve is practically linear. The strain 
data from figure 7 were used to calculate the deflections, y ( - , from deflection equation (11) for the case of 
n = 8. Table 4 compares the deflections, y \ , calculated from deflection equation (1 1) with those obtained 
from the SPAR displacement output. Again, the difference between the values calculated from the SPAR 
program and those calculated from equation (11) is minimal. 


Table 4. Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (1 1); four-ply composite cantilever tube; n = 8. 


Deflection ( y ; - ), in. 



FO 

FI 

F2 

F3 

SPAR 

0.0000 

0.0749 

0.2577 

0.5421 

Deflection 
equation (11) 

0.0000 

0.0658 

0.2519 

0.5413 

Difference , 
percent 

0.0000 

12.2264 

2.2545 

0.1531 


F4 

F5 

F6 

F7 

F8 

0.9110 

1.3475 

1.8347 

2.3558 

2.8946 

0.9167 

1.3610 

1.8570 

2.3877 

2.9356 

0.6301 

1.0000 

1.2176 

1.3526 

1.4154 


The data from table 4 are plotted in figure 8 for visual comparison. The two deflection curves are 
quite close, indicating the accuracy of deflection equation (11). The deflection curves calculated from 
deflection equation (11) for the four-ply composite tube lie slightly above those obtained from the SPAR 
output. 


Helios Composite Cantilever Tube 

Figure 9 shows the bending strain curve, £/ , calculated from the SPAR output, for the strain-sensing 
stations on the Helios composite tube. In this case, the SPAR strain curve is nearly linear. The strain data 
from figure 9 were used to calculate the deflections, , from deflection equation (1 1) for the case of n = 
8. Table 5 compares the deflections, y,-, calculated from deflection equation (11) with those obtained 
from the SPAR displacement output. Again, the difference between the beam-tip displacement values 
calculated from the SPAR program and those calculated from deflection equation (1 1) is minimal. 
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Table 5. Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (1 1); Helios composite cantilever tube; n = 8. 


Deflection ( y ; - ), in. 



TO 

>1 

T2 

T3 

T4 

T5 

T6 

yi 

T8 

SPAR 

0.0000 

0.0240 

0.0673 

0.1324 

0.2167 

0.3166 

0.4284 

0.5479 

0.6721 

Deflection 
equation (11) 

0.0000 

0.0152 

0.0577 

0.1235 

0.2089 

0.3099 

0.4227 

0.5432 

0.6676 

Difference , 
percent 

0.0000 

36.7126 

14.2263 

6.6853 

3.5818 

2.1381 

1.3307 

0.8669 

0.6725 


The data from table 5 are plotted in figure 10 for visual comparison. The two deflection curves are 
quite close, indicating the accuracy of deflection equation (11). The deflection curves calculated from 
deflection equation (11) for the Helios composite tube lie slightly below those obtained from the SPAR 
output. 


DISCUSSION 

In-flight predictions of the deformed shape of the ultralightweight flying wing require the installation 
of multiple bending and distortion strain sensors, equally spaced on the surface of the tubular wing spar, 
for measuring the bending and distortion strains of the wings. The strain sensor data can then be input to 
theoretical displacement equations for calculating slopes, deflections, and twist angles of the wing spar at 
all strain-sensing stations, including the sensing station on the spar tip. 

Because the accuracy of the displacement (bending and torsion) equations has been validated by 
finite-element analysis, the equations are ready for application. To solidify confidence, however, the 
accuracy of the equations must be checked extensively by testing several types of continuous and 
composite cantilever tubular beams under bending, torsion, and combined bending and torsion. The 
validation experiments are currently under way. 

When strain sensors are installed on the wing spar (or wing box), the displacement equations together 
with strain sensor data can be used to construct the in-flight deformed shape of an aircraft such as the 
Helios flying wing. Because they are lightweight, fiber optic sensors may be the optimal choice for the 
embedded strain sensors on the wing spar surface. The use of a conventional strain gage sensing system is 
impractical, because it has numerous lead wires that add weight, thus it is too heavy for the lightweight 
flying wing to carry. 

The displacement theory that has been developed can also be applied to calculate the deformed shape 
of simple beams and plates (ref. 5). The displacement theories can be extended to cases of nonuniform 
beams, such as tapered cantilever wing boxes and aircraft fuselages (varying cross sections) under 
combined bending and torsion (ref. 5). 

CONCLUDING REMARKS 

Displacement equations were developed for a cantilever tubular wing spar under bending, torsion, and 
combined bending and torsion loading. These displacement equations were expressed in terms of strains 
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measured at multiple sensing stations equally spaced on the surface of the wing spar. The principal results 
are as follows: 


1. The displacement equations developed for the uniform cantilever beam were successfully validated 
for accuracy by classical beam theory and finite-element analysis. 

2. The deflections predicted from the deflection equations for the three types of cantilever tubular 
spars 

agreed quite well with the deflections calculated from the finite-element analysis. 

3. The deflection equations and associated strain-sensing system (such as fiber optic sensors) form a 
powerful tool for in-flight calculations of slopes, deflections, and cross-sectional twist angles of the 
wing spar at any strain-sensing station. These calculated displacements can be used to construct the 
deformed shape of the long-span flying wing. Ultimately, the calculated deformation data can be 
visually displayed in real time for the ground-based pilot to monitor the in-flight deformed shape of 
unmanned aerial vehicles. 
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FIGURES 



EDO 1-0209-2 

Figure 1. Solar powered Helios prototype flying wing during a test flight over the Pacific Ocean. 



Figure 2. Cantilever tubular spar with equally spaced bending and torsion strain sensors. 
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Figure 3. Instrumented cantilever tubular spar under combined bending and torsion. 



/ = 100.5 in. 

Nodes 

3673 

a = 4 in. (Al, four-ply walls) 

Four-node elements 

3636 

t = 0.02296 in. (Al, four-ply walls) 
a = 4.07749 in. (Helios spar caps) 

Three-node elements 

36 

T = 0.17784 in. (Helios spar caps) 



c = a + 0.5 1 




060101 


Figure 4. Helios tubular spar finite-element model subjected to a vertical tip load of P = 100 lb. 
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Figure 5. Bending strains, £„ at different sensing stations, x,, calculated from the SPAR program; 
aluminum cantilever tube; n = 8; P = 100 lb. 



Figure 6. Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (1 1); aluminum cantilever tube; n = 8; P = 100 lb. 
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Figure 7. Bending strains, £„ at different sensing stations, x h calculated from the SPAR program; four-ply 
composite cantilever tube; n = 8; P = 100 lb. 



Figure 8. Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (1 1); four-ply composite cantilever tube; n = 8; P = 100 lb. 
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Strain-sensing station identification number, i 


Figure 9. Bending strains, £„ at different sensing stations, x,, calculated from the SPAR program; Helios 
composite cantilever tube; n = 8; P = 100 lb. 



060119 


Figure 10. Comparison of deflections calculated from SPAR with those calculated from deflection 
equation (1 1); Helios composite cantilever tube; n = 8; P = 100 lb. 
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